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ABSTRACT 

Near ultraviolet observations of WASP- 12b have revealed an early ingress compared to the op- 
tical transit lightcurve. This has been interpreted as due to the presence of a magnetospheric 
bow shock which forms when the relative velocity of the planetary and stellar material is 
supersonic. We aim to reproduce this observed early ingress by modelling the stellar wind 
(or coronal plasma) in order to derive the speed and density of the material at the planetary 
orbital radius. From this we determine the orientation of the shock and the density of com- 
pressed plasma behind it. With this model for the density structure surrounding the planet we 
perform Monte Carlo radiation transfer simulations of the near UV transits of WASP- 12b with 
and without a bow shock. We find that we can reproduce the transit lightcurves with a wide 
range of plasma temperatures, shock geometries and optical depths. Our results support the 
hypothesis that a bow shock could explain the observed early ingress. 
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1 INTRODUCTION 

The detection of hot-Jupiters, which are found to be orbiting ex- 
tremely close to their host star with periods of a few days, is 
contin uing to challenge ou r theories of planet formation and evo- 
lution jWatson et al.ll201 11) . Transit surveys such as Super WASP 
provide great i nsight into the properties of extra solar planets and 
their host star jSeager & Mallen-OrnelasI |2003|) . As the number 
of detected planets has increased, work has started on charac- 
terising their atmosphere s iMad husudhan et al. 201T|; ICroll et al] 
l201ll : lHellingetalj|201 l]) including clo ud coverage and tempera- 
ture structure jVidal-Madiar et alj|201lh . 

Following on fro m the identification the hot- Jupiter WASP- 
12b jHebb et al.|[2009l) . recent HST observations indicate an early 
ingress of its near ultraviolet transit compared to the optical data 
( lFossatietaI]|2010l) . There have been several attempts to explain 
such an early ingress by assuming the presence of absorbin g ma- 
terial around the planet jLai et allboidlVidotto et al.ll201(3l) . One 
possible explanation for the presence of an asymmetry is that heav- 
ily irradiated gas giants such as WASP-12b can fill and even over- 
flow the plari etary Roche lobe iOu et alj|2003l ; llbgui et alj|2010l ; 
iLi et al.ll2010l) resulting in mass transfer from the planet to the star 
jLai et al.ll201ol) . 

One further explanation for this observed asymmetry is the 
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presence of a bow shock. Ividotto et alj ( 1201 Ol) have shown that if 
the relative motion between the planetary and stellar coronal ma- 
terial is supersonic, then a bow shock surrounding the planet's 
magnetosphere could form. Such a shock could compress the lo- 
cal plasma to the densities required to reproduce the early ingress 
observed in the near UV lightcurves. 

The possibility of detecting and characterising planetary 
bow shocks provides a new method to study exoplanetary mag- 
netic fields. If the magnetic field of the star can be deter- 
mined through, f or example, Zeeman-Doppler imaging techniques 
jPonati & Landst reel 2009), then the early ingress of the near ultra- 
violet trans it can be used to pla ce limits on the exoplanetary mag- 
netic field ( IVidottoetaLll201lh . The magnetic field is believed to 
provide a key role in shielding the planetary atmosphere from en- 
ergetic particles. 

In this Letter, we present theoretical light curves for a planet 
surrounded by a magnetospheric bow shock transiting its host star. 
Our aim is to test the hypothesis that the early ingress of the near 
UV transit lightcurve can be explained by the presence of such a 
bow shock. 

2 THE MODEL 

If the relative velocity of the planetary and the stellar coronal mate- 
rial is supersonic then a bow shock could form (see Figure[Tll. The 
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Figure 1. Schematic of the shock geometry as viewed along the stellar ro- 
tation axis (not drawn to scale). The shock normal makes an angle ifo to the 
direction of motion of the planet. The distance from the planet to the shock, 
tm, is determined primarily by the strength of the planetary magnetic field. 
Xm denotes the maximum distance between the planet and the projected 
lateral extent of the shock. This distance depends on the chosen values for 
tm, ifio and Ai^. Arm is the region of compressed stellar material behind the 
shock. 



interaction between planetary material and stellar material com- 
presses the local plasma to produce a region of higher density 
plasma behind the shock. If the optical depth of the shocked ma- 
terial is high enough then starlight will be absorbed and produce an 
early ingress in the transit lightcurve. 

The angle between the shock normal and the orbital direction 
of the planet is given by <fo. The value of tpo is determined by the 
relative velocity of the plane tary and stellar coronal material. Fig- 
ure 1 o f lVidottoetai] ( l201(]|) illustrates the scenarios leading to the 
various shock orientations. There are two limiting cases: an "ahead- 
shock" (^0 0) forms when the planet is embedded in the stellar 
corona and a "dayside-shock" (<fo — > 90°) forms when the radial 
wind velocity is very much greater than the relative azimuthal ve- 
locity of the planet. 

Here we use models for the stellar corona and wind 
( IVidottoetai]|20ld. l20T 11) to obtain a value for the plasma den- 
sity at the planet. These models as sume a typical solar base density 
of no ~ 10^ cm"^ l lWithbro3l 19881) and either an isothermal hydro- 
static corona or an isothermal thermally driven wind. We assume 
an adiabatic shock with a maximum compression ratio of 4. For 
an isothermal corona of temperature T, the density of stellar mate- 
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For the stellar wind case we use mass conservation {nii,r^ - const) 
and the momentum equation. 
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to obtain values for the radial velocity, u,. and density, n^,,. The 
plasma density can then be converted into a density of fully ionized 
magnesium using the relation. 



Table 1. Parameters used in our simulations. In all cases we use Xm = 
5.5Rp. The columns are respectively: Model number; Temperature of stel- 
lar plasma; Angle of the shock normal; Calculated value for the number 
density of Mgll; Required thickness of shocked material from our simula- 
tions to reproduce the early ingress; Angular extent of the shock; Calculated 
maximum optical depth for each model. 
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where, nMg/nn is the ratio of Magnesium number density to Hy- 
drogen n umber density whi ch is derived from the metallicity of the 
host star jHebb et alj|2009l) . For WASP-12, MMg/nn = 6.76 x 10"^ 
jVidotto et ai] |2010|) . From this density we can then find bow 
shock geometries a nd orientations that fit the HST observations of 
iFossatietaLlbOld) . 

T o investigate whether the model presented bv IVidotto et al.l 
( l2010h is able to reproduce the data from the near-UV observations, 
we use Monte Carlo radiative transfer calculations to produce simu- 
lated light curves. The parameters we adopt to match the WASP-12 
system are: Mp = \AlMj, R,, = i.l9Rj (where Mj and Rj are the 
mass and radius of Jupiter), M^. = 1.35Mo and R^ = l.STRg. The 
host star is a late F type and the planet orbits in the equat orial plane 
with an impact parameter b = 0.36/?* ( iHebb et al]|2009h . 

The shocked material is considered to be at a distance r^j from 
the planet, with a thickness Ar^ and an angular extent 2Aip. The 
projected lateral extent of the shock is dependent on r/n, <fo and Aif. 
The maximum distance between the planet and the projected lateral 
extent of the shock, Xm, can take the following forms 
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2.1 Monte Carlo Radiation Transfer 

Our simulated transit light cur ves are produced using a 3D Monte 
Carlo radiation transfer code jWood & Revnolds|[l999l) . The cir- 
cumplanetary density structure is prescribed on a 3D spherical po- 
lar grid (coordinates r, d, ip) and is externally irradiated with Monte 
Carlo photon packets with a distribution that reproduces the spatial 
intensity distribution of a limb-darkened star. We assume a spher- 
ical planet and a limb-darkening law such that the intensity / is 
given by 
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where // = cos (? = (1 - r')''-, < r < 1 is the radial distance into 
the stellar disk normalized to the stellar radius and /(O) is the emer- 
gent intensity at the centre of t he star ( Man del & Agol 20021). The 
coefficients a„ are chosen from lClarei] l2004h to match the H-band 
limb-darkening of the host star. We assume that the material ab- 
sorbs or scatters radiation out of the line-of-sight with no scattering 
into the line-of-sight, which is valid for the optical depths required 
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Figure 2. Lightcui-ves and mid-transit images for our simulations. The top panel shows the modelled transits (optical - green; UV - black) with the HST 
observations in red. The black lines represent the results of our simulations. The bottom panel is the mid transit image for each of our models 1A-2B (from 
left to right respectively). 



to produce the early-ingress transits. For this letter we assume the 
bow shock is of uniform density and that the material is static, how- 
ever, our models are very general and can incorporate any density 
structure: analytic, tabulated, or from dynamical simulations. 

2.2 Analysis of the HST observations 

Our goal is to determine the range of shock geometries that can 
provide both an early-ingress and sufficient optical d epth in Mgll 
to acc ount for the absorption seen in the HST data l lFossati et al.l 
bold) . The HST observations (shown in the top panel of Figure 
O enable us to constrain the allowed shock geometries. The first 
data point at phase = 0.86 lies on the continuum of the lightcurve 
meaning the early ingress is yet to begin. This allows us to find a 
maximal value for Xm- If this value is too large, the dip in emission 
in our simulated lightcurves will begin too early. We find that we 
require Xm < 14Rp to ensure a fit to this data point. 

Similarly, the second data point at phase = 0.92 provides us 
with a minimum value for Xm- Since this data point no longer lies 
on the continuum, the bow shock must have started transiting the 
star. The minimal value for the projected shock distance to fit this 
data point is Xm > 4.2i ?„ i Lai et al.l201(]|) . We note that the scenario 
discussed bv lLaietal.li201(]t) where material is accreting from the 
planet onto the star would result in an accretion stream transiting 
before the planet and could therefore also produce an early egress. 

For a given density, the depth of the transit lightcurve is deter- 
mined by the area of shocked material as projected onto the plane 
of the sky. The data point at phase = 0.98 therefore allows us to 



constrain Aru and Atp. If the projected area is too large then the 
lightcurve will be too deep and if it is too small then the lightcurve 
will be too shallow. 

If there is shocked material still transiting after the planet has 
moved olf the stellar disk then we will see a late egress in the tran- 
sit lightcurve. The data point at phase = 1.05 coincides with the 
optical transit suggesting this is not the case. We therefore require 
tpo + A(p < 90° to ensure the near UV transit ends at the same time 
as the optical transit. 



3 RESULTS 

We find that an acceptable fit can be achieved for many differ- 
ent plasma temperatures. These temperatures determine the sound 
speed and the wind speed. They therefo re fix the values of and 
also the density of the shocked material jVidotto et al.l2O10h . 

For certain plasma temperatures, however, we are unable to 
provide a fit to all the data points simultaneously with reasonable 
shock parameters. For T = 1 x 10* K, the calculated density is too 
low and the lightcurve is too shallow. Similarly, for T = 3.93x lO^K 
the density is too high and the lightcurve is too deep. We therefore 
choose to concentrate on two temperatures between these values, 
T = 2 X lO^K and T = 2.5 x lO^K. Because the magnetic field 
geometry is very unconstrained, for each of these temperatures we 
present a solution where the planet is embedded in the corona and 
hence the shock is an "ahead shock" with tpQ = 0° and also one 
where it is immersed in the stellar wind meaning (fo is dependent on 
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the plasma temperature. As illustrative examples we choose hif = 
80° for models lA and 2A and = 40° for models IB and 2B as 
these values will ensure the end of the near UV transit will coincide 
with the end of the optical lightcurve whilst still providing us with 
a large projected shock area. For all cases we then vary Ar^ to find 
a fit to the observations. 

Table[T]shows the parameters adopted for our illustrative cases 
where a fit to all the data points has been found. Figure[2]shows the 
simulated lightcurves and mid-transit images for our models. From 
the simulations it is clear that a range of diiferent shock geometries 
and orientations are able to provide a fit to the data suggesting there 
is some degeneracy in the solutions. 

Firstly we find that Xm is a degenerate quantity. The value 
of Xm is determined by the values of rM,<po and Aip and therefore 
different combinations of these values can produce the same lateral 
projected shock extent projected on the plane of the sky. We have 
found that X^ = 5.5Rp provides an acceptable fit to the data. 

As the temperature of the stellar plasma increases, the wind 
and static coronal models predict a larger density at the planetary 
orbital radius. Therefore the projected area of the shock must de- 
crease to ensure the transit is not too deep and can still fit the data. 
This area is determined by the angular extent A(fi and the radial 
extent Ar^. For the cases where Aip = 80° (models lA and 2A), 
a smaller value of ArM is required to fit the data compared to the 
cases where A(fi = 40° (models IB and 2B). This is a consequence 
of the line-of-sight distance through the shock increasing as (fig de- 
creases, allowing more starlight to be absorbed at the shock front. 

In the optical transit the system is completely symmetrical 
about the mid-transit point, however the addition of a bow shock 
breaks this symmetry in the near UV transit. This can be seen in the 
simulated lightcurves as they are not centred around phase = 1 but 
offset by an amount proportional to Xm- Again if better time sam- 
pled observations could be taken, this offset could be used to pro- 
vide further insight into the stand-off distance between the shock 
and the planet. 



4 DISCUSSION 

Our simulation s have shown that it i s possible to reproduce the HST 
observations of lFossati et alj | |2010|) by assuming the presence of a 
bow shock around the planetary magnetosphere. Using models for 
the stellar corona and wind (with a solar base density) we have 
calc ulated the density of Mgll in the shocked material. 

iLai et al.l ( 1201 Oh calculated the column density of Mgll around 
WASP-12b to be > 1.4x lO'^cm"^. For this calculation they assume 
an optical depth t = 1 in the absorption line of Mgll, a velocity 
V X 100kms~ ' and a characteristic length scale 5 = SRp. From this 
IVidotto et alj | |2010|) found the required number density of Mgll to 
be nMgii <: 400cm"' to reproduce the observed early ingress. Us- 
ing these assumed values we have calculated the extinction cross- 
section of Mgll to be: 

o-uai = = 6.5 X lO-'-'cml (6) 

The maximum optical depth t^.^^ can be then be calculated 
using Equation [6] by setting 5 to be the maximum path length in 
the line-of-sight through the shocked material along with the cor- 
responding number density of magnesium, nMeii (from TablelTJ for 
each of our models. These values are shown in the final column of 
Tabled 



We have found that it is possible to reproduce the observations 
with both lower and higher Mgll densities and that the resultant 
shocked material does not need to be optically thick. 

If similar bow shock structures could be observed in other ex- 
oplanetary systems, transit observations could be useful to p robe 
the presence of planetary magnetic fields. l\^dotto et al] ( l201ll) pro- 
posed a list of candidates that could provide signatures of an early- 
ingress, based on the list of available transiting systems in Septem- 
ber 2010. Should UV observations be obtained for these candidates 
we could apply the model developed here to constrain the allowed 
geometries and orientations for bow shocks. 
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